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Abstract. The aim of this paper is to introduce a class of Hamil- 
tonian autonomous systems in dimension 4 which are completely 
integrable and their dynamics is described in all details. They have 
an equilibrium point which is stable in some rare elements of the 
class and unstable in most cases. Anyhow, the eigenvalues of the 
linearization at the equilibrium point are imaginary and no motion 
is asymptotic in the past, namely no solution has the equilibrium 
as limit point as time goes to minus infinity. In the unstable cases, 
there is a sequence of initial data which converges to the origin 
whose corresponding solutions are unbounded and the motion is 
slow. So instability is quite weak and perhaps no such explicit ex- 
amples of instability are known in the literature. The stable cases 
are also interesting since the level sets of the 2 first integrals inde- 
pendent and in involution keep being non-compact and stability is 
related to the isochronous periodicity of all orbits and the existence 
of a further first integral. 



1. Introduction 



Let us have a quick look of our main results, the details and the 
proofs are in the paper. We introduce the Hamiltonian system in R 4 



(1.1) 



qi =P2 

<?2 = Pl 

pi = -g'{<ii)<i2 
yp2 = -g(qi) 

where g is a C 1 function near and satisfies g(0) = and g'(0) > 0. 
The Hamiltonian function for the system is 

(1.2) H(q,p) = p x p 2 + g(q 1 )q 2 

and another first integral is 

f<?i 



(1.3) 



K(q,p) 



Pi 



v( qi ), V( qi ) 



g(s) ds. 
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The gi,£>2-plane is invariant and the origin is a center on it, namely 
an open neighbourhood C of (0,0) in this plane is the union of peri- 
odic orbits which enclose (0,0). We restrict the phase space to M — 
{(Qi,Q2,Pi,P2) G r4 : (Q11P2) G C , (q 2 , Pi) G M 2 }, so we have a global 
center in the qi,p2-plajie. The Hamiltonian vector fields QVH(q,p) in 



(1.1) and nVK(q,p) are both complete on M (with f2 skew-symmetric 



matrix, see (4.4)). The functions H,K are in involution. The vectors 
V-ff (q, p) , VK(q, p) are linearly independent at each point of the set 
N := {(q,p) G M : (^1,^2) 7^ (0,0)} which is invariant for both the 
Hamiltonian vector fields QVH and QVif. 

A nonempty component r C N of a level set of if, If is diffeomorphic 
to S 1 x 1 and there are coordinates <fi mod 2n and z on S 1 x ffi. such 
that the differential equations defined by the vector field QVH on F 
take the form 

(1.4) <fi = lo , z — v (u,v = const) . 

If (xq, 0) belongs to the projection of T on the gi,p2-plane and T(xo) 
is the period of the first component of any integral curve of QVH on 
T, then v — if and only if T'(^o) = 0. 

A necessary and sufficient condition for the stability of the origin of 



(1.1 ) is that all periodic orbits in a neighbourhood of (0, 0) in the qi,p2- 
plane have the same period so the center is (locally) isochronous. In 
this case all orbits of (1.1) in IR 4 , with (^1,^2) near (0,0), are periodic 
and have the same period. Moreover, a further first integral appears 
as we shall see using a result of Barone and Cesar [3]. This happens 
for instance for the following two functions 

(1.5) g(x) = 1 1 , g(x) = l+x- 



VI + 4x ' w (1 +x) 3 ' 

In Section|2]we show a theorem on isochronous centers taken from 
We reproduce here the proof to be self-contained. Our result relates 
these functions with the even functions and so it shows that isochronic- 



ity in the gi,p2-plane, and Lyapunov stability for (1.1), is a rare phe- 



nomenon among the class (1.1). Isochronous centers were first studied 
by Urabe [5J with a different approach, in [6] a new characterization 
by means of the so called 'involutions' h was found which permits to 
construct the functions giving isochronous centers. 



Generally we have instability of the equilibrium for (1.1) and it is 
very easy to show explicit functions for it, for instance 

(1.6) g(x) = sinx , g{x) = ax + (3x 2 + 72 s , 

with a > 0, (fl, 7) 7^ (0, 0). This kind of instability is quite weak since 



the linearization of the system (1.1 ) at the origin in IR 4 has the double 



eigenvalues ±z^/ ' g' '(0). Moreover, there are no asymptotic motions to 
the equilibrium, namely no motions with the origin as limit point as 
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t — > — oo. In the present paper, instability without asymptotic motions 
is called weak instability. 

In our Hamiltonian systems, instability occurs because there is a se- 
quence of initial data which converges to the origin whose correspond- 



ing solutions are unbounded and the motion is slow, indeed (1.4) shows 
that the coordinate z is an affine function of time. 

Acknowledgement. The picture in Section 3 was made using the appli- 
cation Mathematica by Wolfram Research Inc. by means of the package 
CurvesGraphics6 by Gianluca Gorni of the University of Udine avail- 
able on www . dimi . uniud . it/~gorni. 

2. Isochronous oscillations 
Let us start from the following equation 

(2.1) x = -g(x), g(0) = 0, </(0)>0, 

where g is continuous in a neighbourhood of in R and differentiable 
at (in Section [l] the function g was C 1 , while in the present Section [2] 
the existence of g'(0) is enough) . This o.d.e. has the first integral of 
energy 

(2.2) G(x,x) = — + V(x), V(x)= g(s)ds. 

2 Jo 

By means of this first integral, we easily see that each Cauchy problem 



for (2.1) has a unique solution if g(x) = only at x — 0. 

The potential energy V is a C 1 function and there exists V"(0) = 
g'(0) > . We can restrict the domain of g to an open interval J 3 
such that V is strictly increasing on J n R + , strictly decreasing on 
Jnl_, and for each point x G J there is a unique point h(x) G J with 

(2.3) V(h(x)) = V(x), sgn(/i(x)) = -sgn(x), 

where sgn(a;) is the sign of x. We check at once that the function 

(2.4) u(x) := sgn(x) v / 2^(x) 

is a C 1 diffeomorphism onto the image / = u(J) which is a symmetric 
interval 



«'(0) = V^(0), u G Diff 1 (J;7), 
G J, G /, y G / -y G /. 

Moreover, since h(x) = m _1 (— u(x)) we have that h is also a diffeomor- 
phism: 

heDiS\j;J), h(h(x))=x, 

^ ' ' h(o) = 0, h'(o) = -1. 

We call h the involution associated with V. Remark that the graph of h 
is symmetric with respect to the diagonal which intersects at the origin; 
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indeed (x,h(x)) has (h(x),x) as symmetric point and this coincides 
with the point (h(x), h(h(x))) of the graph. 

Theorem 2.1 (Determining isochronous centers). Let V be a C 1 func- 
tion near ml with V(0) = V'(0) = 0, assume there exists V"(0) > 0, 
let J be an open interval as above and let h G Diff 1 (J; J) be the involu- 
tion associated with V . Then all orbits of x — —V'{x) which intersect 
the J interval of the x-axis in the x, x-plane, are periodic and enclose 
(0, 0). Moreover, they all have the same period if and only if 

(2.7) V(x) = (x-h(x))\ xeJ. 

In this case we say that the origin is an isochronous center for x = 
-V'{x) . 



Proof. By composition with the inverse of the diffeomorphism (2.4), 
the first integral (2.2) gives 2 G[u~ 1 (y), x) = x 2 + y 2 . The first part 
of the thesis follows at once. 

Now, consider a periodic orbit in the x, i-plane which intersects the 
x-axis at xo G J, xq > 0, then h(xo) < is the other intersection. By 
the energy conservation we get the period of the orbit as 

(2.8) 7> ) ' <h ' 



*(*o) yj2(V(x )-V(x)) 
If yo <E I, yo > 0, xo = W 1 ^), then h(xo) = u~ 1 (—y ), and the period 
. u ~ 1{yo) dx 

T(y 

(2.9) 

= 2 



"-M-w) yj2(V(x ) - V(x)) 
yo (u- 1 )' (r) dr f yo ({u- 1 )' (r) + {u' 1 )' (-r)) dr 



wo 



The change of integration variable s = arcsin(r/?/o) gives for yo G I, 
Vo > 0, 

r /2 / \ 

(2.10) T{yo) = 2 J ( K (u- 1 )'(yosms) + (u- i y(-yosms))ds. 
So 

(2.11) lim T(y ) = 2tt (n' 1 )' (0) = 2n/ v / V"(Q) . 

All orbits have the same period if and only if T(yo) equals this limit 
value for all yo G /, yo > 0. By (2.9) this is equivalent to the following 
condition for all z G I, z > 0, 

2 7r _ o f z {(u- 1 )' (r) + (u~')' (~r)) dr 



y/V"(0) Jo Vz 2 - r 2 
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Figure 1. 



We multiply both sides by zj \/y 2 — z 2 and we integrate from to 

y e I, y > o, 

z dz 



71 



(2-13) 



VvW) Jo Vy 2 ^ 2 

* y zdz r z ((u- 1 )' (r) + (u- 1 )' (-r)) dr 
o vV - z 2 Jo Vz 2 - r 2 



By the change of the integration order, see Figure [TJ we get 



(2.14) 



y rv 



o </r 



dz [[u' 1 )' (r) + (u^ 1 )' (-r)J dr 



/y 2 — z 2 \/z 2 — r 



-v — 

= 7r/2 



(2.15) 



2y 



V^(0) 



M Hi/) ~ M X (-l/)> 



yel, y>0. 



We deduced this condition from (2.12), and now we see at once that it 



implies (2.12). By plugging y = u(x) into (2.15) we have 
2u(x) 



(2.16) 



x — h(x), 



x £ J, x > 0. 



This condition is equivalent to (2.7) 



□ 



Formula (2.7) corresponds to formula (6.2) in [6], the proof is in- 
cluded in the proof of Proposition 1 in [6] as a particular case. In [6], 
Section 6, there is the construction we are going to show and also an- 
other one. 

From the theorem we easily get the following result. 

Corollary 2.2 (Constructing isochronous centers). Let h : J — >■ J be 

a C 1 function on an open interval JCR containing which satisfies 
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the conditions in (2.6). Let u> > and define 

,2 



(2.17) V{x) = y (x - h{x)) A , xeJ. 

Then there exists V"(0) = u 2 and all orbits of x = —V'(x) which 
intersect the J interval of the x-axis in the x, x-plane, are periodic and 
have the same period 2tt/u. To get h as above, we can just consider 
an arbitrary even C l function on a (symmetric) open interval which 
vanishes at 0, then a it/ 4 clockwise rotation of its graph gives the curve 
y = h(x) provided the interval is small enough to get a graph. 

For instance, starting from x i— > x 2 / v2 we can (even explicitly) 
calculate 

(2.18) h(x) = l+x - y/l + 4x , V(x) = y (-1 + Vl + 4x 

Notice that the original quadratic function must be restricted to a suit- 
able interval to get a function after rotation. Another construction of 
the potential energies of the isochronous centers is given in [U] (Section 
6 from formula (6.3)). Finally, a simple example is also 

(2.19) h{x) = , V(x) = ^ x 2 i 1 ' 



1 + x v ' 8 \l+x 

Now, let us see some necessary conditions to have constant period 
which can be used to check whether a given function V, or its derivative 
g, locally gives an isochronous center or not. In the sequel V has as 
many derivatives as necessary. We saw how h is related to an even 
function. The even functions have vanishing odd derivatives at 0, so it 
will not be a surprise to see that the derivatives of the involution h are 
not arbitrary. Consider the relation h(h(x)) = x, perform 2 derivations 
and calculate at taking into account h(0) =0 and h'(0) = — 1. In this 
way we get an identity which shows that h"(0) can take any value. By 
3 derivations and calculating at we get h^ 3 \0) = — 3h"(0) 2 /2. Going 
further we see that the even derivatives are free while the odd ones are 
uniquely determined by the preceding derivatives. The first 5 terms in 
Taylor's formula are given by 

(2.20) h(x) = -x + ax 2 - a 2 x 3 + bx 4 + (2a 4 - 3ab) x 5 + o(x 5 ) 



where a, b are free parameters. Next, formula (2.7) shows that the 
derivatives of the isochronous V at are constrained to obey some 
conditions. By means of (2.20) we can calculate V^(0) and V^(0) by 



the lower order derivatives and get 2 necessary conditions. Of course 
we can go forward to infinite conditions. 

Corollary 2.3 (Necessary conditions). Let V admit V^(0) and sat- 
isfy V(0) = V'(0) = 0, V"(0) > 0. Moreover, let the origin be an 
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isochronous center for x = —V'(x) . Then 

(2 21) V W(0) ~ 5y(3)( ° )2 V^(0) - ^ (3) (0)V- (5) (0) UOV^jOr 

To illustrate the necessary conditions let us use it to prove the well 
known lack of isochronicity of the simple pendulum 'x + sin x = 

5sin( 2 )(0) 2 5V^(0) 2 



(2.22) V (4 \0)=sm 
Another example is 

(2.23) V(x 



(3) 



(0) = -M0 



3sin'(0) 3V"(0) 



a 



7 



(a > 0) 

for which the second condition in (2.21) implies /3 = 0, and the first 
condition with (3 = gives 7 = too. 



We could also prove the above formulas in (2.21 ) using the approach 



of Barone, Cesar and Gorni [2] where the first 2 derivatives of the period 
function T{xq) at are computed with a procedure which carries on 
to higher-order derivatives. With some regularity on V they find the 
following formulas for V"(0) = 1 

(2.24) T'(0) = , T"(0) = - A ( |f (3) (0) 2 



4) (O)j . 



3. The dynamics in the Lagrangian framework 

In this section we deal with the following 4-dimensional system de- 
fined by the Lagrangian function L 

L(x,y,x,y)=xy-g(x)y, 

x= -g(x), y = -g'(x)y, g(0) = 0, g'(0) > 0, 

where g G C 1 near in R. This system of differential equations has 
two first integrals G(x, x) and F(x, y, x, y) 



(3.2) 



G(x, x) 



i: 2 



V(x) 



V(x) 



9{s 



ds. 



F(x,y,x,y) = yx + g(x)y. 



The first differential equation (3.1) separates and its dynamics was 



studied in the previous section. We restrict our attention to an open 



interval J as in Theorem 2.1 and to the orbits in the x,i-plane which 
intersect J; their union is an open neighbourhood C of (0,0). Let us 
fix x G J, x > 0, and denote by t X(t, x ) the periodic solution 
of x = —g{x) with (xq, 0) as initial condition at time 0. Next, we plug 
X(t, xq) into the second differential equation in (3.1 ) and get the linear 
equation with periodic coefficient, Hill's equation, 



(3.3) 



V = -9 (X(t,x )) y. 
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The partial derivatives of X(t, xq) give 2 independent solutions of (3.3) 



as one see by derivation of the first equation for X(t, xq), in particular 

d 2 dX , . d d 2 X , 

— — (t,x ) = -= —(t,x ) = 

/„ a\ ot 2 dx dx at 2 

d dX 
= — {-g {X{t, x ))) = -g' {X{t, x )) — (t, x ). 

OXq OXq 

Let us define 

dX 1 dX 

so that 

0(0) = 1, 0(0) = 0, v(o) = o, ^(0) = 1, 

#)0(t)-0(t)v(t) = l. 

Let us denote by T(x ), briefly r, the period of t i-> X(t,x ), then 



t (—7- 0(t + r) is also a solution of (3.3), so a linear combination of 
and ijj, and the previous values at give 

(3.7) <j>(t + r) = 0(r) 0(t) + 0(r) = + <j>(r) ij>{t) 



where we used 0(r) = 1 which comes from the last equality in (3.6) 
with t = r if we notice that V ; ( r ) = and V( r ) = 1- Taking the 
derivative, calculating at t = r, and taking into account ip{nr) = 1, we 
have 

(3.8) 0(t+r) = 0(t)+0(r)^(t) 0((n+l)r) = 0(nr) + 0(r), 
for all fiGZ. This last result is equivalent to 

(3.9) 0(nr) = n0(r), Vn 6 Z. 

If 0(r) = the solution is periodic, otherwise it is unbounded. A 
necessary and sufficient condition to have the former case is 0(r) = 
namely 

(3.10) -gL(T(, ),, ) = 0. 

To go ahed we need the differentiability of the period function 

Proposition 3.1. The period T{xq) of t h-> X(t, Xq) is a C l function 
on {x G J : x > 0}. 



Proof. By formula (2.10), for x E J, Xq > 0, we have 

(3.ii) r(z ) = 

tt/2 



= 2y ^(m 1 )' (u(xo) sin s) + [u l )' (— u(xq) sin s)j ds. 

The function u is a C 1 diffeomorphism as in the previous section. In the 
stronger hypothesis of this section, g 6 C 1 , we have that u(x) = y/V(x) 
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Figure 2. 



for x > is C 2 as well as V = g' and its inverse u 1 (y) for y > 0. This 
proves the result. □ 

Now, we can differentiate d t X(T(xo), x ) = with respect to x 
(3.12) 4^ (T(x ), x ) T(x ) + S^L (T(x ),x ) = 0. 



dt 2 ' dx dt 

We just remind that (d 2 X/dt 2 ) (T(x ),x ) = —g{x ) < for x > 0, 



and have that condition (3.10) is equivalent to 
(3.13) r'(x ) = 0. 

This is the condition in order <fi to be periodic. Since the linear combi- 
nation of (f> and ip gives the general solution, we have just proved the 
following 

Proposition 3.2. Let us fix Xq G J with xq > 0, then the solutions to 



the T{xq) -periodic equation (3.2) are either all T{xq) -periodic, or all 
unbounded but those proportional to t t-> dtX(t,xo). A necessary and 
sufficient condition for the former case is (3.13). 



In Figure [2] the curve (x(t),y(t)) = (X(t,x ),(j)(t)), namely the so- 
lution to system of differential equations (3.1) with x(0) = xq > 0, 
y(0) = 1, x(0) = y(0) = 0, in a non-periodic case. The initial point 
(xo, 1) and the point (h(xo) , g(xo) / g(h(xo))) are crossed at any period. 

Now, notice that the origin in IR 4 is an unstable equilibrium for (3.1 ) 



if and only if we can find Xq > arbitrarily close to such that equation 



(3.2) has unbounded solutions. So by Proposition 3.2 we have 



Theorem 3.3 (Stability and weak instability). Let g be a C 1 function 
near in 1 with g(0) = 0, g'(0) > 0, then the origin in M 4 



stable equilibrium for the system (3.1), if and only if the origin in 
is a locally isochronous center for 'x = —g(x), in this case all orbits of 



is a 
2 
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(3.1) with (x,x) near (0,0) are periodic and have the same period. If 
the equilibrium is unstable, then there is a sequence of initial data which 
converges to the origin whose corresponding solutions are unbounded. 
The instability is weak, by this we mean that there are no asymptotic 
motions to the equilibrium, namely no motions which tends to the origin 
as t — >■ — oo. 

We already proved the first part of the statement, for the last, let 
us only remark that the distance of an orbit from the origin in R 4 is 
greater than the distance of the projection in the x,x-plane which is 
strictly positive unless the projection is (0, 0) and in this case the second 



differential equation in (3.1) gives the harmonic oscillator y = —g'(0)y. 



An explicit function g to give such an example is sin(x) (see (2.22)). 
However, as we know from Section |2j almost all functions g as in the 
statement of the theorem, give instability. 

Examples of the rare functions which give stable equilibria are (see 
pl8] ) and pl9| )) 



(3-14) = \ (l - ^ ) , 9 '(0) = 1, 



1 

4 V" ' ~ (1 + z) 
For all these functions which give stable equilibria, at least whenever 



(3.15) g ( x ) = -(i + x - 7 —^) ) g'(0) = l. 



g G C 2 the system (3.1) admits a (positive definite) Lyapunov function 
which is a further first integral; this additional first integral is smooth 
in a neighbourhood of the origin in R 4 and at least continuous at the 
origin as proved by Barone and Cesar in [3] . In the following statement 
C C R 2 is the set defined below formula ( |3.2[ ). 

Theorem 3.4 (Lyapunov functions). Let g : J — > R be a C 2 function 
on the open interval JC1 with G J, g(0) = 0, g'(0) > 0. Suppose 
that all orbits of x — —g(x) which intersect the J interval of the x-axis 



in the x,x-plane, are periodic and have the same period 2tt/ yg'(0), so 
the origin in R 4 is a stable equilibrium for the system (3.1 ). Then there 
exists 

(3.16) E(x,y,x,y) = a(x,x)y 2 +b(x,x)yy+c(x,x)y 2 , a, b, c : C — > R, 
continuous and positive definite function, which is a (global) first inte- 



gral for the system (3.1). 



The first integral is obtained by means of the functions in (3.5) and 
a suitable inverse function, see [3] for details. 
Let us also notice that 

Proposition 3.5 (Eigenvalues). Let g be a C l function near in R 
with g(0) = 0, g'{0) > 0, then the linearization of the system (3.1) at 
the origin in R 4 has the double eigenvalues ±i^/g'(0). 
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Let us remark that we can easily construct unstable cases with some 
xq at which ( |3.13 ) holds and so the corresponding orbits are all periodic. 
If we want the period function composed with u^ 1 to be for instance 
T(y ) = 2%(1 — Uq + t/q), whose derivative vanishes at y = l/\/2 we 
integrate 



(3.17) 



.</<> 



zT(z) 



So we get the relation 



--dz = 2tt 



(3.18) 



u 



u 



-yo) 



yo 



- 3 vl + Vh vl 



The symmetric choice u 1 {- 
(3.19) vr\y Q ) 



y ) = -u 1 (y ) gives 

= yo - 3% + Y^ y °- 



By inversion we get u{x), then V(x) = u(x) 2 /2 and finally g(x) = 
V'(x). There is an isolate x at which (3.13) holds. We can also 
imagine more complicate examples where these Xq accumulate at 0. 

Finally, let us mention that most of the results in Section [3] were 
found in [7] by a different (more complicated) proof which included 
other differential equations in a family for which the present system 
was a particular case. However, in [7j I had not realized the Lagrangian 
character of the present case which is very important to our purposes 
and will be exploited in the next section. 



4. Complete integrability 
In this section we write q = (q\, q^) = (x, y). The Legendre transfor- 



mation takes the Lagrangian system (3.1) into 

H(q,p) = pip 2 + g(qi)q2, 

dH , , dH 



(4.1) 



dpi 

Pi = -—(q,p) = -g'(qi) 92, 



?2 = -7^— (q,p) =pi, 

op 2 



P2 



dH 



(q,p) = -g{qi)- 



oq\ dq 2 

The Hamiltonian function corresponds to the first integral F while the 
first integral G becomes 

'■'-11 



(4.2) 



K{q,p) 



P\ 



+ V( qi ), V(qi 



g(s) ds. 



Our phase space is 

(4.3) M = {( qi ,q 2 ,p u p 2 ) ER 4 : (q u p 2 ) E C ,(q 2 , Pl ) ER 2 } 

where C C R 2 is an open set mentioned in Section [3] between formulas 
(3.2) and (3.3) so to have a global center in the qi,p 2 -plane. The 
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Hamiltonian vector fields 
/0 



(4.4) QVH(q,p) 







-1 



(4.5) QVK(q,p) 



/o 





1 







-1 







o\ 
1 





V 







o\ 
1 





(d qi H{q,p)\ 
d q2 H(q,p) 
d Pl H(q,p) 

\d p2 H(q,p)j 

fd qi K{q,p)\ 
d g2 K(q,p) 
d pl K(q,p) 
\d P2 K(q,p)J 



t 



P2 
Pi 

\ -9(Qi) I 



( 



\ 





P2 

■g{qi) 
o 



\ 



are both complete on M, namely their integral curves are all defined 
on the whole E. Indeed, this is clear for (4.4) since we have a global 
center on the qi,p 2 -p\ane and linear equations in q2,Pi', for (4.5) we 
simply integrate and remind (4.3) 



(4.6) 



<7i(o), g 2 (t) = g 2 (0)+p 2 (0)t, 
Pl (t) = Pl (0) - g( qi (0)) t, p 2 (t) = p 2 (0). 
The functions H, K are in involution, indeed their Poisson brackets 



vanish: 



(4.7) 



{H,K} 



dH 
dq 

g(<n) 



OK 
dp 


P2 



dK dH 
dq dp 

'gili 








The vectors VH(q,p),'VK(q,p) are linearly independent at each 
point of the set 

(4.8) N:={(q,p)eM:( qi ,p 2 )^ (0,0)} 

which is invariant for both the Hamiltonian vector fields QVH and 
QVK. Indeed, the condition aSJH(q,p) + (3VK(q,p) = implies 
a fl'(?i) — and ap2 = 0, so a = since (gi,p 2 ) ^ (0,0); then it 
also implies (3g(qi) = and (3p 2 = so (3 = 0. 

Let r C iV be a nonempty component of a level set of (H, K), then 
we see at once that it is not compact since it contains the unbounded 
curve pi !->■ (gi, c/g(qi),pi, 0) where c is the value of H and V(qi) ^ 
is the one of K. By a well known theorem (see Theorem 3, Chapter 4, 
in Arnold, Kozlov and Neishtadt |lj), V is diffeomorphic to S 1 x R. 
More precisely 

Theorem 4.1 (Complete integrability) . Let g be a C 1 function near 
in E with g(0) = 0, g'(0) > 0. The functions H,K : M -»■ E, in plj ) 
and (4.2), are in involution on M, dimM = 4, and the Hamiltonian 
vector fields QVH and QWK are complete on M. The set N in (4.8) 
is invariant for QVH and QVK and H, K are independent on N . Let 
V d N be a nonempty component of a level set of (H, K) then it is 
diffeomorphic to S 1 x E and there are coordinates <fi mod 2tt and z on 
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S 1 x K such that the differential equations defined by the vector field 
QVH on T take the form 

(4.9) (p = co , z = v (cj, v = const) . 

If (xq, 0) belongs to the projection ofT on the q\,p2-plane and T(xq) is 
the period of the first component of any integral curve of QVH on T, 
then uo = 2ti/T{xq) and v = if and only ifT'{xo) = 0. 

We already proved everything in the statement of the theorem but 
the last sentence which is a straight consequence of T'(^o) = being 
the necessary and sufficient condition for the solutions on T to be all 
periodic (see Proposition 3.2), and (4.9) gives periodic solutions if and 
only if v — 0. 

For the isochronous systems, namely whenever the period function 
is constant and all integral curves of QVH in M have the same period, 
we know from Theorem 3.4 that a further first integral exist which is 
quadratic in (72, Pi- Of course we do not expect this first integral to be in 
involution with K. In the next section we investigate the (isochronous) 
cases where a third first integral quadratic in the momenta pi-,P2 exits 
and we arrive at explicit formulas for it. 

5. Explicit superintegrable systems 

Let g be a C 1 function near in M with g(0) = 0, g'(0) > 0. Suppose 
that the system (4.1) has a first integral W of the form 

(5.1) A(q u q 2 )p\ + B(q u q 2 )piP2 + C(q u q 2 )pj + U (qi,q 2 ) ■ 

Then the equation {H, W} = give a cubic polynomial in pi,p 2 - The 
coefficients vanish if and only if 

(5.2) d 2 A = 0, diC = 0, 8^ + 0^ = 0, d 1 B + d 2 C = 1 
diU(q u q 2 ) = 2C(q u q 2 )g(qi) + B(q u ^M^Ofe 
d 2 U(q ll q 2 ) = 2A(q l ,q 2 )g'(qi)q 2 + B(q 1 ,q 2 )g(q 1 ). 



(5.3) 



The conditions in (5.2) are the same as those in (3.2.2) in Hietarinta jl] 
where the celebrated Darboux problem is treated (only the order of the 
coordinates is changed since gi = p 2 , q 2 = p\ ) . They have the following 
solution where a, b±, b 2 , cj, c 2 , C3 G R are arbitrary 

^(91,92) = aql + b x qi + c x 
(5.4) B{q u q 2 ) = -2aq x q 2 - b x q 2 - b 2 q x + c 3 

C(QuQ2) = aqj + b 2 q 2 + c 2 



Plugging (5.4) in (5.3) we get the following condition for the integra- 
bility of U 

(5.5) 3(b 2 + 2aq 2 )g( qi ) = 

= 3(6i + 2aq 1 )q 2 g'(q 1 ) + 2{c x + g 1 (6 1 + aqx))q 2 g" {q x ) . 
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For q 2 = we get 3b 2 g(q 1 ) 
we have 



which implies b 2 = 0, and back to (5.5) 



(5.6) 
For a 

(5.7) 



0, 9(0) 



3(&i + 2aq 1 )g'{q 1 ) + 2( Cl + q^h + aq 1 ))g"{q 1 



and g'(0) = u: 

,2 



to 

T 



1 this equation gives 
1 



1 



there are no solutions) 



The values 



where A = &i/(2ci) (for ci 
A = 2 and a; = 1 give (3.14). The first integral (5.1) for this g is easily 
obtained with b\ = 2Aci, a = 0, b 2 = 

(5.8) cip?(l + 2Xqx) + pip 2 (c 3 - 2\ Cl q 2 ) + c 2 pj+ 

+ c 2 ^ (-1 + a/1 + 2Agi) + Ciw 2 g2 + ?2 (c 3 - 2Acig 2 ) g{qi) + d 

where del. For c\ = C3 = d = 0, c 2 = 1/2, we have the first integral 
K we already know, while H corresponds to the choice c\ = c 2 = d — 0, 
C3 = 1. For C\ = C2 = 1, C3 = d = we have the first integral 



(5.9) 



P?(l 



2Agi) - 2Ag 2 p lP2 +p\ + 2V( qi ) + (w 2 - 2\g{ qi )) q\ . 



This function is positive definite near the origin of 1R 4 as we check at 
once. This is why we have Lyapunov stability of the equilibrium and 
compact orbits. 

We dealt with the particular case a = 0. In the sequel a ^ and 
we fix a = 1 dividing the first integral by a. The coefficient of g" in 
equation (5.6) is a quadratic function in qi with discriminant b\ — Ac\. 
Let us consider first the case where this vanishes, namely c\ = b\/A. If 
b\ = then (5.6) with g(0) = 0, g'(0) = oo 2 , gives the trivial function 
g(x) = u 2 x, while for A := 2b± ^ we get the following solution with 



2. 



g(Q) = 0, g'(0) = u 
(5.10) g( gi 



to 

T 



A + gi 



A J 



(A + gi) 3 

Notice that this formula includes (3.15) for u — 1 and A 
integral (5.1), which corresponds to (5.10), with C3 
which vanishes at the origin of M 4 , is 

(5.11) p\(X + qi f - 2p lP2 (X + q x )q 2 + p\(\ + g 2 2 ) + 



1. The first 
0, c 2 = 1, and 



+ 



(x + qi f + x* 1 + 4q2 



X 2 u 2 



(A + gi) 2 . 

We can check that at the origin of 1R 4 its gradient vanishes and the 
Hessian matrix is the diagonal matrix 2(u 2 , X 2 u 2 , A 2 , 1), so we have 
again a (positive definite near the origin) Lyapunov function for all 
values of the parameter A 7^ 0. 
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Now we are ready to deal with the generic case. Let a = 1 and 
b\ — 4ci 7^ 0, we also consider b\ ^ and c x ^ since b± = gives the 
trivial solution and C\ = no solutions at all. Then the solution of the 



equation (5.6) with g(0) = and g'(0) = or 2 > is 



(5.12) g( qi ) = {b 2 C ^ Ci)2 1 {b\ + 4c 1 )(b 1 + 2 gl ) + 



fc 6?-4c 1 -2(6 1 + 2g 1 ) 2 ) 
+ Oi / 

1 + ?i(6i + <?i)/ci 



The first integral (5.1) with c 3 = 0, c 2 = ci, divided by c 1; and which 



vanishes at the origin of Mr, is 

(5.13) pi +pl + — (pigi -V2q-i) {pi{h\ + gi) -^2^2) + 

Cl 

2 



+ 



26iW 2 (&! + 2 9l )(-4ci( Cl + (61 + q x )q x ) + (b 2 - A Cl )q 2 



- 4c i) 2 A/l + gi^i + gO/ci 

We can check that at the origin of 1R 4 its gradient vanishes and the Hes- 
sian matrix is the diagonal matrix 2(u 2 ,u 2 , 1, 1), so we have a positive 
definite function near the origin. 
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